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Recent numerical simulations of a disordered system J3J/ have shown the existence of 
two different relaxational processes (called stimulated and spontaneous) characterizing the 
relaxation observed in structural glasses. The existence of these two processes has been 
claimed to be at the roots of the intermittency phenomenon observed in recent experiments. 
Here we consider a generic system put in contact with a bath at temperature T and char- 
acterized by an adiabatic slow relaxation (i.e. by a negligible net heat flow from the system 
to the bath) in the aging state. We focus on a simplified scenario (termed as partial equi- 
libration) characterized by the fact that T = (where only the spontaneous process is 
observable) and whose microscopic stochastic dynamics is ergodic when constrained to the 
constant energy surface. Three different effective temperatures can be defined: a) from the 
fluctuation- dissipation theorem (FDT), T e g DR ; b) from a fluctuation theorem describing 
the statistical distribution of heat exchange events between system and bath, T e g T and c) 
from a set of observable- dependent microcanonical relations in the aging state, T^§ R . In a 
partial equilibration scenario we show how all three temperatures coincide reinforcing the 
idea that a statistical (rather than thermometric) definition of a non- equilibrium temper- 
ature is physically meaningful in aging systems. These results are explicitly checked in a 
simple model system. 

1. Introduction 

Efforts to extend well established concepts in equilibrium thermodynamics to the non- 
equilibrium domain have repeatedly appeared many times in the past in different contexts. 
An idea that has attracted the attention of physicists for quite a long time is the concept 
of a temperature applied to non-equilibrium states pQ. In equilibrium, the notion of 
temperature can be covered from two different perspectives. On the one hand, there is the 
thermodynamic approach where temperature is defined as a parameter that characterizes 
classes of systems in mutual thermal equilibrium. The usefulness of this thermometric 
definition relies on the validity of the zeroth law of thermodynamics. On the other hand, 
there is a statistical approach (ensemble theory) where temperature can be defined from 
the properties of individual systems without any reference to the mutual equilibrium 
property. The statistical temperature is defined as the inverse of the energy gradient 
of the configurational entropy measured over a constant energy surface. The maximum 
entropy postulate relates the statistical notion of the temperature to the thermometric 
one. The statistical and thermodynamic concepts look equivalent but they are not as the 
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latter requires a specific behavior of the different systems when put in mutual contact. 
The thermodynamic definition of a temperature represents a stronger condition than the 
statistical one. 

An interesting category of systems that has recently received considerable attention 
are glassy systems in their aging regime. The aging regime is a slow relaxational process 
characterized by the extremely small net amount of heat delivered by the system to the 
bath in contact. It has been suggested [2] that in aging systems a thermo metric definition 
of non-equilibrium temperature is meaningful for a thermometer responding to low-enough 
frequencies. We suspect that this definition is too strong and might be wise to investigate 
a low-level statistical definition of a non-equilibrium temperature. 

The purpose of this paper is to show that a statistical definition of non-equilibrium or 
effective temperature (rather than thermometric) can be rescued for a particular category 
of aging systems characterized by adiabatic slow relaxation. In this category of systems 
relaxation is imposed to be ergodic when constrained to a given energy shell. At T = 
this condition defines what we term as partial equilibration scenario. For this class of 
systems three different definitions of an effective temperature are possible: a) from the 
fluctuation-dissipation theorem (FDT), T^ DR ; b) from the fluctuation theorem applied 
to the statistical distribution of heat exchange events between system and bath, T e g T ; c) 
from a set of microcanonical relations in the aging state, T e ^ R . In a partial equilibration 
scenario we show how all these temperatures coincide reinforcing the idea that a statistical 
definition of a non-equilibrium temperature is physically meaningful in aging systems. 
We check all statements and results in a simple model of glassy system where the partial 
equilibration scenario holds and explicit computations can be done. 

2. Adiabatic relaxation 

Let us consider a system that is prepared in a non-equilibrium state by placing it 
in contact with a thermal bath at low temperature (quenching protocol). The system 
will relax and release heat to the bath in a process that can span from several minutes 
to millions of years. Only when the system has reached equilibrium the net heat current 
from the system to the bath vanishes. All along the paper we will refer to this relaxational 
regime as the aging regime and the corresponding non-equilibrium state as the aging state. 
The time after the quench will be referred as the age of the system and will be denoted by 
one or two of the variables s and t depending on whether one-time or two-time quantities 
quantities are considered. We adopt the convention t > s. 

Several definitions and quantities seem appropriate to put the discussion in perspec- 
tive 1 . Let us consider a system described by an energy function E{C) where C denotes a 
generic configuration. The system is in contact with a thermal bath and the microscopic 
stochastic dynamics is both ergodic and satisfies the detailed balance property. PtiC) will 
denote the probability for the system to be in the configuration C at time t. The average 
value of an observable A at time t will be denoted by A(t) and is given by, 

A(t)=J2A(C)P t (C) . (1) 
c 

1 For a detailed presentation of several of these concepts see [3] 
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Often the time argument in A(t) will be dropped off and we will simply write A, with the 
clear understanding that in general it is a time- varying quantity. P(C, t\C, s) denotes the 
conditional or transition probability for the system to be at configuration C at time t if it 
was at configuration C at time s. The autocorrelation (C^t, s)) and response (Ra{P, s)) 
functions are defined by, 

C A (t,s) = Y / P(C',t\C,s)A(C)A(C') ; R A (t,s) = (2) 

Both Ca and R A can be decomposed into a stationary (fast) and aging (slow) parts. In 
the large s regime, for many aging systems a quasi-FDT relates the aging parts of (J2J) in 
terms of the fluctuation-dissipation ratio (FDR), 

• (3) 

T e g DR (s) is a time-dependent effective temperature that has been shown to have some of 
the properties of a thermodynamic temperature [5J. In a weak ergodicity breaking scenario 
CU(£, s) decays to zero for t — s — > oo in a typical time that we denote by Tdecorr(s) 2 . In 
general, Tdecorr(s) oc s a with a a given exponent. In many cases a = 1 appears to be a 
very good approximation, usually termed as simple aging. 

We define the aging regime as adiabatic if the fraction of heat released from the system 
to the bath goes asymptotically to zero for time differences of order of the decorrelation 
time, 



E{s + 

Tdccorr 

(s))-E(s) 



E(s) 



for s large enough (4) 



Of course this relation is meaningful only in the aging regime where the energy E(s) is 
still far from its equilibrium value 3 . The relation between the rate decay of the energy 
and the frequency domain where FDT violations are observed has been already pointed 
out in jl]. 



3. Stimulated and spontaneous relaxation as the origin of intermittency. 

For an adiabatic regime as described in (jlj) two type of heat exchange processes are 
predicted to be observable depending on the timescale of the measurement. According 
to (HI), and for times of the order of Td 0C orr(s) or smaller, the net heat flow transferred 
from the system to the bath is exceedingly small, yet heat fluctuations can be as big 
as if the system were in equilibrium at the bath temperature. In such case there is 
a continuous heat exchange between the system and the bath and energy fluctuations 
< E 2 (s) > — < E(s) > 2 (measured over r dccorr (s)) are of the same order of the energy 
content E(s) and determined by the heat capacity of the system at the bath temperature. 

2 This could be defined either as the integral correlation time f CA(t',s)dt' or the value of t — s for 
which CA{t, s) decays to 1/e of its maximum value Ca(s, s). 

3 Deviations from the power law divergence Tdocon (s) oc s a may provide an explicit check whether the 
system has not left the asymptotic regime and aging is not interrupted (in this last case the decorrelation 
time starts to saturate to its value at equilibrium). 
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Were one to measure at age s the probability distribution P S {Q) of heat exchanged Q 
between the system and the bath along intervals of a given duration r ~ T decorr (s), a 
Gaussian distribution would be found with zero mean (as no net heat is transferred from 
the system to the bath) and a variance c 2 (T) which is independent of the age s but 
dependent on the temperature of the bath. This feature is ubiquitous in glassy matter 
(e.g. structural glasses quenched at low-enough temperatures) where the net heat flow 
from the system to the bath is unobservable, yet heat is quickly exchanged with the bath 
as thermal conductivity is high 4 . For all practical purposes the system looks equilibrated 
at the temperature of the bath and a thermometer put in contact with it would measure 
the bath temperature. We call this heat exchange process stimulated as it originates in 
the existence of physical processes thermally excited by the bath. 

For times much larger than r,i ecorr (s) the average net heat flow from the system to 
the bath is not negligible, a clear consequence that the system has not yet equilibrated. 
Most of the exchange processes that occur along these timescales are stimulated by the 
bath, however other exchange events do not fall into the previous category and follow a 
completely different distribution. Contrarily to the stimulated process, this distribution 
is expected to display an exponential tail whose width A(s) depends on the age of the 
system s. However this process is not excluded to be described by a Gaussian distribution 
as well. Indeed an exponential tail is characteristic of Gaussian distributions centered 
around a finite value. We call this new process spontaneous as its statistical description 
is not determined by the temperature of the bath but rather by the fact that the system 
has been prepared in a non-equilibrium state. Figure [T] illustrates typical distributions 
for the stimulated and spontaneous processes. The existence of two different kinds of 
heat exchange process (and therefore two different heat exchange distributions) occurring 
along widely separated timescales has been recently numerically verified in the context 
of a simple spin-glass model [3]. It has been suggested that the existence of these two 
processes is at the roots of the intermittency phenomenon observed in glasses [6j and 
colloids [7j 5 . 

4. Statistical (microcanonical) description of the aging state 

Along the rest of the paper we will analyze in detail and idealized particular case of the 
more general previous scenario. In doing this we aim to understand fundamental issues 
behind a possible statistical, rather thermodynamical, description of the glassy state. We 
will consider in detail a situation where the stimulated process does not exist and only 
the spontaneous process is observed in the adiabatic relaxation. We will refer to it as 
the partial equilibration scenario. It is then possible to prove that effective temperatures, 
with a precise statistical meaning, do emerge. This is the content of the next Subsections, 
I4.1|4.2|4~31 In Sec. 03 we will focus our attention on the oscillator model (OM) |HJ as a 
simple case where all results of the current Section can be explicitly verified. 



4 The simplest evidence in favor of this statement is a piece of silica vitrified at room temperature (i.e. 
not yet equilibrated) whose temperature can be felt by touching it with the hand. 

5 For the latter it might seem more appropriate to speak about stress release rather than heat exchanged 
as compaction (rather than heat dissipation) is the main relaxational process. 



5 




Figure 1. The origin of intermittency: Two mechanisms of relaxation, stimulated (continuous 
line) and spontaneous (dashed line) characterize the overall relaxation. Q stands for the heat 
transferred from the bath to the system (the sign indicating the direction of flow, if Q < 
heat flows from the system to the bath) along regularly spaced time intervals and P S (Q) the 
corresponding probability distribution at age s. The black points illustrate what experimental 
measured data would look like. In the simplest scenario both processes are described by a 
Gaussian distribution. The stimulated component has a variance <r 2 (T) dependent on the bath 
temperature but independent on the age. The spontaneous component shows an exponential tail 
in the Q < side (shown as a continuous line) whose slope A is age dependent. The spontaneous 
process for Q > is unobservable as relaxational dynamics is constrained by the existence of a 
net heat transfer from the system to the bath. 



4.1. Partial equilibration scenario 

To suppress the stimulated process the bath is put at zero temperature (therefore, no 
heat can flow from the bath to the system). However, to keep the system in an adiabatic 
relaxational regime (jlj) (avoiding the aging regime to be quickly interrupted due to the 
presence of forever-lived metastable configurations) we require an additional condition: 
Dynamics must be ergodic when the system is constrained to move in any constant energy 
surface 6 . This condition is illustrated in Fig. El In this case, it is easy to prove that 
relaxation cannot arrest at T = as it is always possible for the system to decrease 
its energy by moving along the constant energy surface until a favorable downhill move 
occurs. Models that fall into this category have been termed as models with purely 
entropic barriers, the OM |5] being a prominent example. Eq. (J3J states that in an 
adiabatic regime the energy remains practically constant during timescales of the order 
of the decorrelation time Td eC orr(s). Therefore, as correlations decay in timescales much 
shorter than the energy does, the system can explore a large number of configurations 



6 A more precise statement requires dynamics to be ergodic within the energy shell between E and E+AE 
where AE is finite in the large volume limit. 
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Ergodic trajectory 
Non-ergodic trajectory 



Figure 2. Schematic picture of the type of ergodic dynamics considered in a partial equilibration 
scenario. Any phase-space points A, B contained in the surface of energy E must be connected 
by a path lying in that surface (dashed line). Note that other paths (continuous line going 
through C) that we call non-ergodic (as they change the energy of the trajectory) might also 
exist as the original dynamics of the system coupled to the bath is ergodic. 



within a given region of phase space. In case only spontaneous relaxation takes place, the 
system is constrained to reach a quasi-stationary state where all configurations lying in a 
given constant energy surface have the same probability. Equiprobability is guaranteed by 
the fact that microscopic dynamics over the constant energy surface is ergodic and satisfies 
detailed balance, i.e. it is reversible over the constant energy surface. The process where 
a given region of the constant energy surface is sampled according to the microcanonical 
ensemble (i.e. all configurations lying in that surface have the same probability to be 
visited) will be termed partial equilibration. 

For later purposes we define the complexity Se{A) by the relation 

S E (A) = \og(n E (A)) = \og[j2 S(E - E{C))8{A - A(C))} (5) 

c 

where Qe{A) stands for the number of configurations with energy E and observable A 7 . 
For the total complexity S(E) we have, 

S(E) =log(n(J5)) =log[J2S(E-E(C))} (6) 

c 

where now Q(E) is the number of configurations with energy E irrespective of the value 
of the observable A, Q(E) = J cIAQe{A). In (jHJ) we use the energy label E as a subindex 
and not as an additional argument, just to stress its key role in the dynamics compared 
to that of other observables. Note that we could equally well use the subindex s, as for a 
given age s the average energy of the system is fully determined by the dynamics. 

7 As we are counting configurations, Sg(A) is nothing more but the usual configurational entropy. How- 
ever, we deliberately avoid to use this term and prefer to talk about complexity. We have in mind the 
most general case (not addressed in this paper, see instead 0) where partial equilibration is established 
among regions or components of phase space. For this case the use of the term configurational entropy 
could be misleading. 
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4.2. A fluctuation theorem in the aging state 

In a partial equilibration scenario transitions between configurations lying at different 
constant energy surfaces (i.e. transitions that increase or decrease the energy) are not 
equiprobable. However, after applying an external perturbation, a shift of energy levels 
takes place and some configurations, initially belonging to different energy surfaces, may 
end up into the same one. We show below in Sec. 14.31 how the response of the system 
to an external perturbation is determined by the density of energy levels just around the 
reference value E(s). Accordingly, transition rates between configurations having different 
energies E, E', were they equiprobable, are described by the following microcanonical 
relation, 

P(AE) _Q(E') 



P(-AE) n(E) 

where AE = E' — E is the (intensive) energy difference between both configurations. 
It is important to point out that P(AE) is not the actual transition rate, but the rate 
the system would display if transitions between configurations having different energies 
were equiprobable, a situation that is encountered only after shifting the energy levels (by 
applying an external perturbation) and redistributing them into a unique energy surface. 

Eq. (jZJ) has the form of a FT since it describes the ratio of rates in the forward E —>■ E' 
and the reverse E' —>■ E directions. Fluctuation theorems, similar to (J7J), have been 
derived in other contexts, for instance in systems in steady states jU] or systems arbitrarily 
perturbed from an initial equilibrium state jTUj. From © we can write, 



P(AE) \(dS{E) 



AE 



p(-ae) ~ ex n(^) AE n = ex nf^ii)\ (8) 

where we have expanded the complexity S(E) and kept only the first term in AE 8 . The 
factor T e g T (s) in the exponent in the r.h.s of (fTTj) defines an effective temperature, 

1 f S(E), 



T e F ff T (s) V dE )e=e { s) (9) 

where we specifically include the subindex s to denote its time dependence through the 
value of the energy E at age s. The super-index FT indicates that this effective tem- 
perature is derived from the fluctuation theorem (JJJ). For a Gaussian P(AE) the value 
of 

^eff T ( s ) can De shown to be proportional to the width of the exponential tail A(s) as 
depicted in Fig. [TJ This connection has been exploited in |Hj as a possible way to estimate 
the effective temperature from intermittency measurements. 

Eqs. (1718)1 may look as simple detailed balance but it is not. Detailed balance is a in- 
herent property of the microscopic dynamics. Two are the main differences between ()8I9|) 
and microscopic detailed balance. In the former we assumed equiprobability between 
transitions with identical energies and this is not guaranteed if not in the partial equili- 
bration scenario. Moreover the factor in the exponent of (jHJ) is not the temperature of 



8 As AE is an intensive quantity, higher order powers in AE should be included in JHJ. However they are 
not relevant for what is addressed in the present paper. A similar situation is later encountered in (J8J. 
See the ensuing footnote © for a more detailed explanation. 
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the bath as implied by detailed balance, but a quantity that is solely determined by the 
dynamics. 

We contend to show that T e g T (s) as derived from the fluctuation theorem coincides 
with the effective temperature derived from the FDT relations that link correlations and 
responses in the aging regime. The origin of this connection has been already mentioned. 
Eq. (j7J) links transitions between configurations with different energies. These transitions 
can be probed only by lifting the energy of the different configurations after applying an 
external field. 

4.3. Microcanonical relations and effective temperatures in the aging state. 

Let us consider A to be any observable of the system that can take different values in 
a given surface of constant energy E. The equiprobability assumption implies that the 
transition rates between configurations with different observable values A, A', at a given 
age s (when the surface explored has energy E), satisfy the following relation, 

W E (AA) Q E (A') 

We(-aa) n E (A) 1 ' 

where A A = A' — A and Q E {A) was defined in ((HJ). The ratio of rates fllOD is therefore 
age dependent as the value of E changes with the age of the system. Eq. (fTU|) says that 
the rate for the observable A to change its value in AA within the surface of constant 
energy E, when going from the value A to the value A', is proportional to the number 
of configurations with value A' at the energy E. Although (fTUj) describes rates there is 
no explicit reference to any timescale. In fact, we use the term W for these rates (as 
compared with the term P(AE) used in 0) to stress the fact that these are rates rather 
than probabilities (i.e. they have dimensions of a frequency). Note that the difference 
between using P or W is minor as any timescale drops off when considering the ratio 
among the forward and reverse rates. Again, as was done for (JJJ), we expand the term 
in the r.h.s. of (fTU|) around A' = A by using (jSJ) and consider only the first term in the 
expansion, 

W E (AA) w dS E (A)- 
- exp ' 



W E (-AA) 



(^fP)AA + 0((AA) 2 )] (11) 



The variation AA = A' — A is intensive, therefore all powers of AA enter in (jllj) even in 
the thermodynamic limit. However, as we will see later, only the first term is relevant for 
the emergence of effective temperatures. Here we will not discuss the possible relevance 
of higher order terms 9 . Relation (jllj) says that dynamics evolves towards configurations 
that have a higher value of the complexity S E (A), i.e. towards the maximum value A*(E) 
oiS E (A), 



>dS E {A) 



= (12) 



V QA >A=A*{E) 

As the value of E is time dependent it is also the value of the maximum A*(E). In a 
partial equilibration scenario, A relaxes (along timescales of order of Tdecorr(s)) to the 

9 Higher order terms in AA are expected to contain information about the stability of the aging state and 
thermally induced fluctuations of the effective temperature. Contrarily to the temperature of the bath 
(which cannot fluctuate) the latter might display fluctuations even in the large V limit. 
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value A*(E) given in (|12j). This value changes in time as the energy decreases and is 
slaved to the time evolution of the energy. As the energy decreases much slower than A 
relaxes to A*(E) (c.f. (J3J)), reversibility holds in the aging regime, 

W E (AA) = We(-AA) . (13) 

An interesting and special class of observables are those called neutral where, for large 
enough times, the value of A*(E) is independent of E (i.e. of the age of the system) 
A*(E) = A*. Therefore the value A* must coincide with the equilibrium value A eq if the 
system has to equilibrate. Neutral observables have the interesting property that they 
reach a stationary value exponentially fast and their time evolution is not slaved to that of 
the energy. The most prominent example of neutral observables is the wave vector density 
Pk(t) in supercooled liquids which stays negligible at all times as there is no long range 
order that develops in the amorphous glass state. The dynamics of a neutral observable 
is therefore quite easy to visualize. Starting from any initial configuration the dynamics 
of the system quickly evolves towards A* and stays there forever. 

Non- neutral observables are expected to relax fast to their value A*(E) and display 
an interesting non-monotonic behavior as the value A*(E) changes with the age s. For 
instance, if the system starts far from equilibrium but with a value of A corresponding 
to the equilibrium value A eq , initially the system will depart from this value, and follow 
the time evolution of A*(E) to eventually come back again to A eq after equilibrating. 
This effect has been observed in structural glasses where the volume density displays 
non-monotonic behavior in the glass state and is known as the Kovacs effect 10 . 



<A(t)> 



non-neutral behavior 



neutral behavior 



Figure 3. Time evolution of neutral and non-neutral observables. Explanation is given in the 
text. 



We are now in a position to understand the emergence of effective temperatures as 
usually derived from fluctuation-dissipation relations. Suppose now that at time s an 
external field of intensity h coupled to the observable A is applied to the system. In the 
presence of a field the energy of a configuration C is shifted by the Zeeman term, 

E{C) -> E{C) - hA{C) . (14) 
10 This effect has been studied in different models of glasses, for instance |11I12| 
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Under the field, the surface of constant energy E does not coincide anymore with that 
at zero field, meaning that configurations with identical energy at zero field (i.e. lying in 
a given surface of energy E) get shifted by different amounts in a field (i.e. finish into 
different surfaces with different values of E). In particular, for a given value of the energy 
E, from all configurations initially lying in the surface at h — 0, after switching the field 
some configurations leave the surface, others come onto the surface, finally others stay 
there (i.e. those with A(C) = 0). 

According to ()14|). just after the field has been switched on, configurations with a larger 
value of A decrease their energy relative to their energy at h = and configurations 
with a lower value of A increase their energy. However, the complexity S(E) © is a 
monotonic function of the energy, therefore configurations with higher energy are more 
numerous than those with lower energy. As a result of the action of the field, the number 
of configurations Q E (A) with a given value of A, that lie inside the surface of constant 
energy E, monotonically increases with A, 

n h E {A) = J dE'{l E ,{A)5{E -E' + hA) = Q E+hA {A) = 

n E (A) [l + 9 J^Sl A h + 0(h% = n E (A) exp[(gMd) + 0(h 2 )} . (15) 
Using the equivalent of ©, S E {A) = log(Q E (A)), we get, 

S h E {A) = S E (A) + (?^l)Ah + 0(h 2 ) = S E (A) + ME, A) Ah + 0(h 2 ) (16) 
where we have defined, 

ME,A) = (^^) (17) 

V 8E ' E=E(s),A=A(s) K ' 

As S E (A) is a monotonic increasing function of E, ()15j) indicates that relaxation in a 
field is pushed towards configurations with progressively increasing values of A. In a field 
partial equilibration occurs again in the new surface of energy E as all configurations 
there contained remain equiprobable. The new (jlOj) reads, 

Wt(AA) n E (A') W E (AA) , n /rn ASAllS , \ 

where fll5|lfjjl have been used. A word of caution in the derivation of (|TH|) is in place. We 
have considered (3 e e(E,A) = f3 e s(E,A'). This is justified as the quantity S E (A) entering 
into the definition (JT7j) is extensive, i.e. proportional to the volume V of the system, while 
A A = A' — A is an intensive quantity. In the large V limit, the partial derivative (JT7j) 
is the same whether it is taken at A or A' as the difference AA = A' — A is intensive so 
p eS (E,A) - (3 eS (E,A f ) ~ 0(1/V). 

Eq. (J18)) is reminiscent of detailed balance, however the same remark has to be made 
here as was done in Sec. 14.21 The quantity /3 e g(E, A) is not the temperature of the bath 
anymore but a time dependent value as E(s) and A(s) change in time. The quantity 
fl e ff(E,A) defines an effective temperature as obtained from the microcanonical relation 

(MR) (HHD, 

1 - ME ,A) = . (19) 



T™ R ' A (s) ' V dE / E=E(s),A=A(s) 
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As A(s) = A*(E(s)) (see (|T2|l and the ensuing discussion) we can replace the complexity 
Se(A) by the total complexity S(E) in the r.h.s of (|19p. Notably this gives a value of 
T e g R (s) that is independent on the type of observable so the label A in the l.h.s of (fTTJJ) 
drops off, 

T MR (s) PeS W { dE ) E=E[S)A=A{S) { dE ) E=E{S) > l^j 

which coincides with the result obtained from the FT (JHJ). 

From relation (fTH|) it is now possible to link the autocorrelation and response func- 
tions (J2J through via the FDR, where the effective temperature fl2Uj) appears explicitly, 

^(M) ^ ( / C f' 3 W -*) • (2D 

We are not going to show here the details of this derivation but only mention the main 
steps. The proof follows the scheme of the derivation shown in Section 3 of Ref. [B] for the 
standard derivation of FDT where ()18p here corresponds to (49) there. In [3| it is shown 
how the response function can be decomposed in two different contributions called 
and R( 2 \ In equilibrium, when FDT holds, vanishes and R^ gives the equilibrium 
response. In the partial equilibration scenario, (J21j) is obtained whenever R^ vanishes. 
Close inspection of (53) in j3] shows that this term vanishes only if the perturbing field does 
not modify, to linear order in h, the trapping time distribution measured over transitions 
starting at the surface of energy E. More precisely, consider a sample of all trapping 
times for configurations of age equal to s. As dynamics is stochastic, the same quenching 
protocol will generate different configurations with values E, A in the vicinity of their 
dynamical averages E(s),A(s). Let p s {r),p^{r) be the trapping time distributions both 
at zero field and after applying a field h at time s respectively. Eq. (|21|) holds whenever 
Ps(t) = Ps(t) + 0{h 2 ) 11 . In other words, up to linear order in h the sole effect of the 
field is to modify the density of configurations as indicated in ()15|) . 

The physical significance of these results can now be appreciated at its full extent. 
In the simplest partial equilibration scenario the effective temperature derived from the 
microcanonical relation (|2(J|) is independent on the type of observable and coincides with 
that obtained from the FT (jHJ) and the FDR Q. For the last equality, the trapping time 
distribution is required to remain unchanged at linear order with the intensity of the field 
h. In this case, 

Tj s DR (s) = Tl s T (s) = TT(s) (22) 

This equality is remarkable as it shows that the effective temperature can be directly 
obtained from the fluctuation theorem Q without the necessity of introducing a per- 
turbing field and measuring correlations and responses 12 . The first equality in (J22|) has 
been numerically verified in a given example of spin-glass model 0. In what follows we 
exemplify these results for a simple solvable case. 



11 Similarly one could say that the average trapping time is not modified to linear order in h 

12 From a different perspective, a result where a perturbing field is not required to measure the FDR has 

been recently proposed |13I14| . 
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5. The oscillator model (OM) 

Exactly solvable oscillator models (OMs) [H], unrealistic as they look, provide a simple 
physical basis to describe glassy dynamics. A parallel can be established between oscillator 
models for glassy dynamics and the original urn models by the Ehrenfest's. The former 
can enlighten our comprehension of the essential features of glassy dynamics, in the same 
way urn models have provided a ground basis to understand key concepts of equilibrium 
thermodynamics such as the Boltzmann entropy. The OM is an example where the partial 
equilibration scenario holds for the dynamics at T = 0. It provides an excellent framework 
to verify the statements of previous sections. 

5.1. Energy relaxation 

The OM |8j is described by a set of continuous variables Xi and an energy function, 

K N 

£ = y£*? (23) 
z i=i 

where K is the spring constant and N is the total number of oscillators. Oscillators 
are non-interacting and therefore the model has trivial static properties, the total energy 
E = NksT '/2 according to the equipartition law. We consider a cooperative dynamics 
where oscillators are updated according to the rule x,- L — > + the rj being uncorrelated 
Gaussian variables with r\ = and variance r\rj = A 2 5ij. The updating of all oscillators 
is carried out in parallel and the moves are accepted according to the Metropolis rule. We 
will focus our analysis on the dynamics at T = where only updates that decrease the 
total energy are accepted. This leads to slow dynamics as most of the proposed moves 
tend to increase the energy and only few of them are accepted (i.e. the acceptance rate is 
quickly decreasing with time). Dynamics in the OM is ergodic if confined to a constant 
energy surface (see Fig. EJ). Therefore, dynamics does not become arrested at T = as 
no metastable configurations (except the ground state) exist at T = 0. 

Dynamical properties in the OM are derived from the distribution of attempted energy 
changes P(AE). There are several ways to compute this distribution [Ej, the simplest 
one derives from the Gaussian character of such distribution. The change in energy of an 
elementary move is given by 

AE = -^=Y x iri + — Yrl (24) 

From the Gaussian character the of r^, it follows that AE has a Gaussian distribution 
whose mean and variance are given by, 

M AE = AE = KA 2 /2, a AE = (AE) 2 - (AE) 2 = 2KA 2 ^ = 2KA 2 e (25) 



which yields jH], 

{AE-^fy 



P{AE) = (AireKA 2 )-* exp 



4eKA 2 



(26) 
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where e = E/N is the energy per oscillator. At T = 0, the dynamical evolution of the 
energy e and the acceptance rate a(e) (the fraction of accepted moves) are given by the 
following closed equations, 



de r° r° 1 KA 2 

— = j ^xP(x)dx ; a(e) = J ^dxP(x) = -Erfc^— ) (27) 

with Erfc(x) = (2/y/Tr) f£°duexp(—u 2 ) the complementary error function. As the energy 
e decreases the variance of the distribution (J23|) decreases. Accordingly, the acceptance 
rate also decreases. The dynamical evolution of the energy and acceptance can be solved 
in the long-time limit, E{t) ~ 1/ log(t), a(t) ~ l/(tlog(t)). 

5.2. Effective temperature 

Correlations and responses have been computed for the magnetization M = Y^i x i |H] ■ 
In the rest of the paper we will consider A = M as the observable under quest. M is a neu- 
tral observable as can be verified by solving the dynamical equation for the magnetization. 
It relaxes exponentially fast to zero which is the equilibrium value of the magnetization at 
all temperatures. The autocorrelation function C^ifi s) = (1/N) J2i x i(t) x i{s), the corre- 
sponding response function R M (t,s) (J2J) and the susceptibility Xm(^s) = J s RM(t,t')dt' , 
do not have stationary part but only aging part and show simple aging with r docorr (s) oc s. 
Correlations and responses define an effective temperature T eS (t,s) through the relation, 

dC M (t,s) 

Tl a DR (t, s) = — (28) 

RM{t, S) 

which is the temperature T of the bath in equilibrium. At T = the response is finite in 
the OM due to the shift of energy levels described in Sec. 14.31 A simple expression can 
be derived in that case, 

T e F r (t, s) = T?™(s) = Ms) + 77^)^ (29) 

where f(s) is a given function jH] which asymptotically decays as 1/s. Two remarkable 
facts emerge from (|2HJ): 1) T c s(t, s) only depends on the lowest time s, therefore charac- 
terizes the aging state of the system at time s; 2) The second term in the r.h.s of (|29|) is 
sub-dominant respect to the first term leading to T e s(s) — > 2e(s), i.e. the equipartition 
law is asymptotically satisfied in the aging regime. 

5.3. The fluctuation theorem 

As the energy decays logarithmically and Td ec0 rr(s) oc s, (JU) is satisfied and relaxation 
is adiabatic. Moreover, dynamics in this model is ergodic if constrained to the constant 
energy surface. Under these conditions, only spontaneous relaxation occurs (Sec. EJ) and 
partial equilibration holds (Sec. 14.1)1 . We can verify the validity of the FT by substitut- 
ing (j2HI) in ©, 

P(AE) /A£ x /A£n . . 

7 5 F ^ = exp( w )=exp(^-) (30) 

giving T e g T = 2e. This result coincides with the asymptotic value T e g DR previously ob- 
tained (EH- 
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5.4. Trapping time distribution 

Here we compute the trapping time distributions without and in a field P s (t),Ps( t )- I* 1 
the presence of a field the energy of the OM is given by, 

K n N 

E = Ne T = — J2x 2 l -hY,x i = N(e-hm) (31) 

where e,m denote (K/2)x 2 and x (i.e. the energy and magnetization per oscillator re- 
spectively), the total energy per oscillator being ey = e — /im. Eq. ((%H) can be written 
as E = yZ)i=i(a ; j — h/K) 2 — Nh 2 / (2K). The updating rule for the shifted variables 
Hi = Xi — h/K remains unchanged. Consequently, the evolution equation for the quantity 
e' = (K/2)y 2 = e - hm + h 2 /{2K) = e T + h 2 /{2K) is identical to that obtained for e 
at h = 0. As the time evolution of e' does not depend on h this implies that ey gets 
corrections that are even powers of h. In a field, (|27|) holds by replacing e with e'. There- 
fore, the acceptance rate is not modified at linear order in h. In general, the trapping 
time probability distribution is given byp^(r) = aVXl-aV))^ 1 wh ere t is the finite 
number of Monte Carlo steps (i.e. does not scale with N) and a h (e') is the acceptance 
rate in a field. The same expression is valid for p s (r) putting h = 0. This immediately 
proves that the distribution of trapping times remains unchanged at linear order in h. In 
particular, the average trapping time is finite and given by the inverse of the acceptance 
rate, 

TO T 



^ I^P h {AE)d{AE) a\e>) (32) 



where r is a microscopic time. 



5.5. Microcanonical rates for the magnetization 

Let us consider the joint probability P h (AE,AM) of having a change in the total 
energy AE and magnetization AM in the presence of an external field h (E is given in 
(|3*Tj) and includes the Zeeman term). At T = 0, the probability that an attempted 
change AM is accepted is given by, 



W^{AM) = —f P h (AE,AM)d(AE) = — P°(AM) [° P^(AE\AM)d(AE) . (33) 



— f P h (AE,AM)d(AE) = r ' 

T J -co T 

where P h (AE,AM) = P?(AE\AM)P°(AM), the subscript c standing for conditional 
probability 13 . Again, it is easy to show that these distributions are all Gaussian. Straight- 
forward calculations give, 

P»(AA/) = (27rA 2 )-iexp[-^Ml] (34) 
P?(AE\AM) = (2 m l)-ie W [- {AE -f 2 - hAM) ) (35) 

with the definitions ah = Ch&h, bh = 9h&li a \ = ^Ke + h 2 — 2Khm, c? = A 2 (l — ( Km ~ h ) \ 
gh = (Km — h)/(2Ke + h 2 — 2Khm). In the linear response regime we obtain, = 



13 



We are following the notation of Sec. 14. 31 where We(AM) stands for the rate at zero field. 
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f e + (m 2 /(2e 2 ) - \/{2Ke))h + 0(h 2 ), c\ = c 2 (l+mh/e + O(h 2 )), c 2 = A 2 (l -m 2 K/{2e)) 1 
<j\ = 2Ke(l - hm/e + 0(h 2 )),b h = mK -h + 0{h 2 ). Inserting flME^ in (J33D gives, 

Wl(AM) = W exp [_^_ (1 + ^)] (36 , 
N AM) = If^Erifcf^^^l (37) 



«0 



where P°(AM) is given by (|34jl . As m is a neutral observable that relaxes fast to we 
can replace m = everywhere in all previous expressions. Up to linear order in h the 
rates are given by, 

W^(AM) = W E (AM) exp{^^) ; W E (AM) = P°(AM)a(e) (38) 

where a(e) is given in ()27p. We remark the following results: 1) the perturbed rates are 
multiplicative 14 , 2) (0 holds as W E {AM) = W E (-AM) and 3) relation (|TEJ) is satisfied 
as well, 



W£(AM) W E (AM) r (AM)h 



W%(-AM) W E (-AM) 



2T^(e) 



(39) 



so the effective temperature here obtained T e ^ R (e) = 2e again coincides with that derived 
from the FDR and the FT ((HDD - 

6. Conclusions 

A statistical interpretation of a non-equilibrium or effective temperature in aging sys- 
tems, rather than a thermometric one, might be possible. Recent numerical simulations 
of a disordered model [5] have shown that the effective temperature, as measured from 
FDT violations, originate from the existence of a spontaneous relaxational process de- 
scribing heat exchange low-frequency events. Superimposed to it there is a stimulated 
process that is characterized by the temperature of the bath and responsible of most of 
the heat exchange observable events between system and bath. Although the timescale 
of the spontaneous process is related to the temperature of the bath, its statistical de- 
scription (e.g. the specific form of the corresponding FT) is related to other properties 
of the aging state such as the energy content. The concurrence of these two processes 
in the overall relaxation is related to the intermittent phenomenon recently observed in 
various experiments |6|7j . Moreover, the width of the exponential tail associated to the 
spontaneous process is predicted to be proportional to the effective temperature. 

We have investigated a partial equilibration scenario where only the spontaneous pro- 
cess occurs (the temperature of the bath is set to zero) and dynamics is ergodic when 
constrained to a given energy surface. In this case three effective temperatures can be 



14 The form of the perturbed rates (|38f) is multiplicative. This assumption was made in a given class of 
trap models at the level of configurations and shown to lead to the existence of effective temperatures [16) . 
Although the multiplicative property may hold at the coarse-grained level of observable values, it was 
shown to be far-fetched at the level of configurations J7| . Eq. i|38[) shows that rates are to be considered 
multiplicative only at the level of observable values. 
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defined: a) from the fluctuation-dissipation ratio Q; b) from the fluctuation theorem for 
the statistical description of the spontaneous component (J7J); c) from microcanonical rela- 
tions relating observable changes (J2(J|) . All three are shown to be identical (|22|). Explicit 
calculations in the OM demonstrate the validity of these statements. 

Several open questions and directions of research along the present ideas appear worth- 
while. It would be useful to go beyond the qualitative level of demonstrations in the 
present paper over more founded mathematical proofs of all these results. Attempts try- 
ing to establish the origin of effective temperatures using master equation formalisms have 
already appeared in the literature |13|18j and we foresee more work in the near future. 
We should also mention the close similarity between the present approach and that by 
Edwards for granular matter [19] . Would be very interesting to investigate other classes 
of models with simple equilibrium properties, such as kinetically constrained models [20] , 
where the existence of non-trivial effective temperatures is still under debate. In these 
models it is possible to identify correlated motion of particles [2T| that underpin a spa- 
tially heterogeneous dynamics [22.- Another category of interesting problems to explore 
are systems in steady states where FDT violations have been studied and where 
non-Gaussian effects, similar to those described here, have been identified as well |24| . 
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